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Abstract 



A simple model is presented for the calculation of the quenched average over 
impurities which are rendered static by setting their mass equal to infinity. 
The path integral formalism of the second quantized theory contains annealed 
averages only. The similarity with the Gaussian quenched potential model is 
discussed. 



I. INTRODUCTION 

Disordered systems can be characterized by two time scales, the time t spent between the 
preparation of the impurity distribution within the sample and the measurement, and the 
impurity diffusion time, t{. There is actually a third characteristic time, the time span of the 
measurement, t m but it is far from the other time scales, t m <<t, U. We are interested in the 
dynamics in the limit when t, ti — > oo but R = t/U stays finite. The initial distribution of the 
impurities within a given sample determined by the preparation remains unchanged when 
R « 1. For R » 1 the impurities reach equilibrium with the faster degrees of freedom, 
electrons. The self averaging quantities measured on large samples can be obtained by 
means of averaging over the impurity distributions. The latter, the impurity distribution is 
determined by the preparation method when R << 1 or by the equilibrium properties for 
R » 1. 

The usual modelisation methods were developed for R « 1. The time available for the 
impurity motion in this case is insufficient to establish the feedback of the fast degrees of 
freedom on the impurity dynamics. Such a suppression of a part of the dynamics can be 
realized by the replica method [§] the introduction of spurious supersymmetric particles 
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or the Keldysh formalism ||J]. A simple model to describe electrons in the presence of 
quenched disorder is based on the partition function 0] 

1 n - J [v]Z5 [■?/; t ]Z) [-0] e~ ^ dx ' d *^i (x,t)(iS t + ^ r A-^Cx))^- (a;,t)] ^ ^ 

where v(x) is a static potential representing the disorder, ipj denotes the electron field, and 
j = 1, • • ■ ,n is the replica index. The limit n — > is made formally at the end of the 
computation of the observables. 

Although these methods have already led to a number of important results they are 
not completely satisfactory. The strategy based on supersymmetry can not cope with the 
Coulomb interactions. The Keldysh formalism, given in terms of the propagators, is inher- 
ently perturbative. The replica method is not always reliable due to its formal nature. 

Apart from glassy materials or fast quenching, the usual preparation of samples gives 
R » 1 and there is time for the electron dynamics to influence the impurity distribution 
seen during the measurement. One can easily incorporate the feedback of the electron 
dynamics into the impurity distribution by avoiding the replica method in the model (ffl), 
i.e. by setting n = 1. But the simple Gaussian distributed static field which is sufficient 
for the modelisation of the forces due to decoupled impurities might be too simple for the 
reproduction of the full impurity dynamics with its non-trivial time dependence. 

We describe in this paper a simple model for R » 1, when the slowly moving impurities 
are in equilibrium with the electrons. Instead of the static potential v (x) and the replica 
method used in (H), we return to the more detailed impurity dynamics. The impurities 
will be made static as far as the observables are concerned, by setting their mass equal to 
infinity. Their equilibrium with the electrons will be established by considering the partition 
function at finite density or temperature. Our model involves particle degrees of freedom, 
the impurities and the electrons. The elimination of the former produces an effective model 
similar to ([!]) with n = 1. We find that the salient feature of quenched disorder, the 
possibility of generating localisation, is realised in the same manner for R » 1 as for 
R << 1. In particular, the perturbative solution of our model for the electronic observables 
obtained by resumming the two-loop self energy can be mapped onto a subset of graphs 
coming from ([]]), namely those graphs which have no more than two impurity lines ending 
at the same space location (our impurities obey fermion statistics). This result implies 
that the conductivity, when computed through the resummation of the maximally crossed 
diagramms in the particle-hole channel, agrees in the two models. Since our model is given 
in terms of an annealed partition function this result opens the way for the application of 
different non-perturbative numerical methods to deal with quenched disorder when R » 1. 

The organization of the paper is the following. Section II introduces our model which 
contains impurities with infinite mass, i.e. zero mobility. The infinite mass non-relativistic 
particles have a singular Fermi sphere, reflecting the high level of degeneracy of the ground 
state. This makes the use of the grand canonical ensemble questionable in this case. The 
canonical ensemble is not spoiled by the degeneracy of the free Hamiltonian. In section 
III we present a method to compute averages in this ensemble. The propagators are com- 
puted in section IV and section V describes the setting up of the perturbation expansion 
in the framework of the path integral formalism. The self energy and the propagator for a 
particle-hole pair are computed in section VII. Finally conclusions are drawn in section VIII. 
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Technical points concerning the regularization and the perturbation expansion are developed 
in the Appendices. 



II. THE MODEL 

Consider the system of electrons and impurities described by the fields ijj and </>, respec- 
tively. The hamiltonian is of the form 



H = I dx 

n 2 



2M 



h 2 A 
4> ] (x)A(j)(x) + gipl(x)ijj a (x)^(x)(j)(x) 



(2) 



where a and b stand for the spin indices. Electrons and impurities are present with finite 
density, what is usually achieved by using the grand canonical ensemble, the modification 
of the Hamiltonian, 



H - dx 



fx e iplip a (x) + Hitf {x)(t>{x)\ . (3) 

Here \x e and /ij are the chemical potential for the electrons and the impurities, respectively. 
The spinless impurity field is supposed to be anticommuting, a reasonable approximation 
even for integer spin atoms in solids. In Eq. (^) H is considered as an effective hamiltonian 
for the conducting band with a spatial resolution which is so low that the inhomogeneity of 
the crystalline structure is not seen. 

We render the impurities static by setting M = oo. Then the lagrangian possesses an 
additional, space dependent U(l) gauge symmetry 

<f>(x,t) ->e*( x >0(z,t), (4) 

which indicates the absence of spatial correlations for the impurities. This is not equivalent to 
the space-dependent part of the electromagnetic U(l) gauge invariance, it only characterizes 
the difference between the slowly moving and the truly quenched impurities. 

What happens in the limit M — > oo? To understand this limit better it is useful to 
introduce the free diffusion constant D = H/2M for the impurities. The impurities generate 
an interaction on the distance scale yDT during the elapse of time T. When this distance 
shrinks well below the spacelike cutoff, the lattice spacing, the model reaches the quenched 
limit as far as the observables up to time T are concerned, and the local symmetry (|J) 
is approximatively realized. The distribution of these apparently static impurities is not 
uniform, it is governed by the ensemble used in the computation. Since we shall use the 
canonical ensemble, our results refer to a thermally equilibrated and after then quenched 
impurity distribution. 

It is easy to see that the limit M — > oo enhances the large momentum and the low 
frequency contributions of the interactions. The perturbative short distance contributions 
are strengthened by the absence of the kinetic energy in the denominator of the propagator. 
This increase of the short range fluctuations is the impact of the static nature of the impu- 
rities on the spatial disorder. The slowing down of the impurity motion induces non-local 
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interactions in the time rather than the space direction ||. The increased sensitivity at low 
frequencies arises from the suppression of the hopping between neighboring spatial locations, 
the reduction of the effective spatial dimensionality of the impurity dynamics <i e // — > as 
M — > oo. This reduces the impurity dynamics to the sum of non- interacting Quantum 
Mechanical single site problems. The result of the strong correlations generated in time is a 
possible non-local effect in time within the electron sector. 

The dependence of the static observables on M/m is not continuous at M/m = oo, 
e.g. the conductivity diverges for M/m < oo because the electrons, moved by the external 
electric field drag the impurities along. The conductivity is finite only in the absence of recoil, 
when M/m = oo. The discontinuity at M/m = oo is reflected in the existence of the gauge 
symmetry (§]) and the singularity of the Fermi surface. In fact, the Fermi sphere for M < oo 
consists of the one particle states p 2 /2M < //j. The Fermi momentum is pp ~ p l ^ d in d 
dimensions where p is the particle density which yields //j ~ p 2 l d /M. Hence the Fermi sphere 
becomes a highly singular point in the momentum space as M — ■> oo and the construction 
of the grand canonical ensemble is not obvious. The same problem is seen by the vanishing 
Fermi velocity, Vp ~ p l l d /M. The ground state is highly degenerate and admits the gauge 
symmetry (f|). Due to this complication we shall carry out the computations below for 
M = oo in the canonical ensemble for the impurities. 



III. CANONICAL ENSEMBLE FOR THE IMPURITIES 

Let us consider first a system of fermionic impurities without spin. In order to avoid 
singularities in the static limit we place them on a lattice with N lattice points (see Appendix 
0) but keep the continuous notation whenever it is not misleading. The hamiltonian with 
an external source J is chosen to be time dependent, 



H (t) 



dx 



h 2 



2M 



4>\x, t)A(f)(x, t) + J(x, t)4>\x, t)(f)(x, t) 



(5) 



in the Heisenberg representation^. We introduce the real time expectation values of the 
operator O in the canonical ensemble corresponding to the density p as 



(0) p = -^—TrPpT 



■ f* T dtBo(t) 



o 



(6) 



where 



TrP p T 



-if* T dw (t) 



(7) 



and P p is the projection operator onto the subspace of n = pV particles (V denotes the 
volume), 



By anticipation of the perturbation expansion in the source J the expressions for the free 
propagator will be given for J = 0. 
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P J-n 27T 



(8) 



which guarantees finite density for arbitrary mass. The finite temperature averages can be 
obtained by performing a Wick rotation in the time evolution, 



(O) 



-TrP n T 



with 



Z PtP = TrP p T 



-J o dtH o (t) o 



(9) 



(10) 



The trace is evaluated by using localized states either in momentum or in real space for 
M < oo or M = oo, respectively, 



TrO = ^2(n\0\n), 

{n} 



(11) 



where the states \n) are the Fock states determined by the occupation number configuration 



n. 



\n) 



n p (at)" P |0) forM<oo, 
U x (ai) nx \0) forM = oo, 



and 



[dp, aj] + — 5 PjQ , [a x , aj+ — S XjV . 



(12) 



(13) 



The products of creation operators are taken in an arbitrarily chosen but fixed order in ([12] 
We obviously have the relations 



(ala q ) p = p5 P:Q , (ala y ) p = pd x . 



(14) 



for a translational invariant system, where p = a d p stands for the dimensionless impurity 
density. These results allow us to obtain the free propagator in the canonical ensemble. The 
causal propagator is defined as (cf. Appendix ||) 



iG(x,t,x',t') = (T[(f)(x,t)(p\x',t')]) 
giving, for a fixed volume V, 



iG(x, t, x', t' 



1 

V 



E 



ip(x— x')— ie p {t— t') 



1 - p for t > t', 
—p for t < t', 



(15) 



(16) 



where e„ is the single particle excitation spectrum for J = and 



iG(x, t, x', t' 



[1 - p)5 X)X > for t > t', 
-f>8 x>x > for t < t', ' 



(17) 



for M < oo and M = oo, respectively. The Fourier transformation 
G(p,q,uj) = 5 Piq G(p,uj) 

roo roc 

iG(p, u) = {l-p) / e^-^dt - p / e - it{uJ ~ ep) dt 
Jo Jo 

is performed with the usual choice of the boundary condition 



e lut dt 



Acut—St 



dt 



(18) 



(19) 



with 5 0+ 



G(p, u) 



G(x, x', uj) 



l-p 



P 



u — e p + i8 uj — tp — i8 



uj — e p + id 



+ 2nip5{uj — € p ) 



uj + i8 

The retarded Green function can be written as 

1 



+ 2nip5(uj). 



Gr(p,uj) 



uj — e„ + i8 



G R (X,X',UJ) = S x>a 



1 



uj + i6 



(20) 



(21) 



and the advanced one is Ga = G* R . 

In order to gain an insight into the impact of the impurities on the dynamics consider 
the coupled impurity-electron system defined by the hamiltonian 



H 



dx 



n 

—il>l(x, t)Ai(j a (x, t) + g${x, t)(f)(x, t)ipl(x, t)ip a (x, t) 
2m 



+-(V>10M)V>aOM) 



(22) 



Since the impurities located at different space locations decouple it is easy to obtain the 
effective theory for the electrons where the impurity degrees of freedom appear as local 
variables. Leaving the actual computation to Appendix [D] the final result for the expectation 
value of an observable O between the electronic states |^/i, e ) and \^ff, e ) i n the canonical 
ensemble for the static impurities can easily be obtained as a summation over the impurity 
occupation numbers n%, J2 x n x —n — pN, N = V/a d , 



(<%e|^,e>> = ^E<%el T 



O 



1^1,6/] 



(23) 



where 



{n-x} 



I * i,e/ ; 



(24) 



and 
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H e ff(t; n$ 



dx 



h 2 

~ q— Vfe t)A^p a (x, t) + gn(x, t)ipl(x, t)ip a (x, t) 



+ T (VlOM^aOM))' 



(25) 



This result shows clearly the connection with model ([]]), namely the impurities are static 
both models but they reach equilibrium with the electrons in our case. 



IV. AUXILIARY PROPAGATORS 



The expressions (^) give the correct propagator in the canonical ensemble but are not 
well suited for a perturbation expansion. In order to generate more conventional Feynman 



rules we look into the detailed way the averages (|T4| ) were achieved and introduce some 
auxiliary quantities. These are obtained by postponing the integration in the projection 
operator (^) and by putting the operator of the exponent in (^]) into the hamiltonian. This 
amounts to the usual strategy of gauge models, treating the integral variables of the con- 
straints in the path integral as dynamical variables. With this in mind we introduce the 
impurity hamiltonian 



H a (t) 



dx 



2M 



f (x, t)A(f>(x, t) + (a+ J(x, t))<p\x, t)(f)(x, t) 



(26) 



with the source J being kept vanishing in this Section and the expectation values given by 



(0) f 



Sl n dae iapv TrT 



O 



JZL dae ia P v TrT 



-| J_ T dtH a ,(t) 



(27) 



where a' = a/2T. It is furthermore useful to introduce the propagator corresponding to a 
given value of a, 



G a (x,t,x',t') = TrT 



One can easily find its explicit form for M = oo, 



G a (x,x',t) =5 x , x ,e- ia,t (l + 



AT-l f 1 



for t > 0, 
-e~ ia for t < 0, 



(28) 



(29) 



since this non-normalized expectation value corresponds to the time evolution of a state 
with a single filled and N — 1 empty sites, see Appendix [B] for the detailed derivation. As a 
check we compute the impurity density, 



-G(x,0,x, 0^ 



J^dae ina (l + e 



,N-l 



(30) 



a d JZ W dae ina (1 + e~ { 

The integration over a selects the particle combinations with the desired particle number 
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N-l 

-G(x,0,x,0 + ) = ^ W ~J/ = S = P- (■>' ) 



When M < oo we find 

1 for t > 0, 

(1 + e >»') <j 

where e p denotes the one-particle energy. The finite temperature propagator reads 



G a (p, t) = n (l + e-*^) { L e — (ep+Q ) ^ J < o' (32) 

si-Art ** 



J^dae<«^<l) a//9 



_ _ ip(x-a: ')-ie p {t-t ') J _^_ **<j7=P V / /oo\ 

7^ dae*™ Ug (1 + e-^- ia ) 1 ; 

which is just the Gibbs average of single particle contributions in the given particle number 
sector. Note that the parameter a is purely imaginary in an imaginary time formalism, it 
appears as the time component of a gauge field. 



V. PERTURBATION EXPANSION 

We have two different realizations of the canonical averages depending on whether the 
projection operator is inserted once or after each dt time stepf] during the time evolution 
— T < t < T. The first case leads to the path integral expression 

Y j da j D^D^Je* f<tefa ^**-^** A * +a ^*-^ +J ***]<9, (34) 

for the canonical average of the operator T[C?], where 

Z p = J da J £)[^]£)[0*] e */*«te[^V-^^A^-a(^^-p)+J^*]. (35) 

In the second case we have a time dependent a(t) trajectory to integrate over, 

— f D [a] D [(f)] D [(f)*] e l J dtdx l i &*&-^0* A <f } + a 'f > * ( i } - a P+ J 'f > * ( i , ]Q_ (36) 

and 

Z P = J Z}[a]D[0]D[0*]e i / d * dx[ ^^-^** A ^ +ar ^- ap+Jr ^. (37) 



2 dt serves as the ultraviolet cutoff needed in the derivation of the path integral formulae. The 
limit dt —* is convergent provided that the number of the degrees of freedom is kept finite by a 
cutoff in space. 



S 



Due to the periodicity of the integrand the integration over a can be performed over the 
whole real axis. The steps leading to fl34j) and ( |36l) are similar to those giving the path 
integral expressions in QED for the static temporal gauge, doA = and for the real temporal 
gauge A = 0, respectively, A (x,t) playing the role of the projection operator parameter 
a(t). 

The straight perturbation expansion for (|36|) yields 



(O) = — [ D[a]D[</>]D[</>*]e i J^te*+-£+*^ a d 

xO^—(i J dtdx{a + J)(jf^j " , (38) 



Z p = J D[a]D[<j>]D[<f>*]e i S^M**-&P A *- a rt 
1 



and 



xj-^i j dtdx(a + J)(j)*(j)j . (39) 

The small parameter of the perturbation expansion is J which stands for the interaction 
with an external source. The small parameter for the projection operator is h for bosons 
when the expansion is performed around <p(x, t) = J~p. 

For fermions there is no simple way of saturating the path integral and the projection 
operator must be implemented nonperturbatively. This can be achieved when it is inserted 
once only in the average, 

(O) = y J da J D[«/,]D[0*] e i /* <fa; ^^-w^ A ^+^^-#P] 

xO Y, -r (i I dtdxJ<p*<f)) " , (40) 



hi 



and 



Z p = J da J D[0] J D[0*] e i / dtda: ^^-^^ A ^>-#^ 



x Y,— ] ( i I dtdxJ<p*<f)) . (41) 



In what follows we shall use this formalism, where the hamiltonian is given by (f26|). 
There is a vertex corresponding to the interaction term J4>*4>, the propagator is G a and the 
integration over a is to be done after the loop integrations. The usual Feynman rules are 
applicable for the computation of the integrand for the integration over a in the numerator 
and the denominator independently. Due to the independent integrations over a the discon- 
nected diagrams do not always simplify in the expectation values, a remnant of the non-local 
nature of the canonical ensemble. But one can verify that the grand canonical result where 
the disconnected contributions simplify is recovered in the thermodynamical limit. 

A short discussion is now in order about the use of the lattice regularization in the com- 
putation of the loop integrals emerging from the perturbation expansion. Non-relativistic 
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quantum field theory is non-renormalizable in itself, a problem which is made even more se- 
rious by sending the impurity mass to infinity as mentioned in Section y. This is naturally a 
formal problem only since the cutoff is actually kept finite in effective theories. Nevertheless 
it is useful to distinguish observables which diverge from those which stay finite when the 
cutoff is removed because the computation of the latter is simpler. To understand the reason 
let us start with the remark that the lattice regulated model is more complicated than the 
one obtained in the continuum due to the trigonometric functions in the propagator. In 
fact, the rule of generating the lattice propagators from the continuum say for electrons in 
2D is 

I > I (42 ) 

u,-& + v + iS& --^(sin 2 ^ + sin 2 ^) +/ , + ^ 2 ' 1 } 

where a s is the spatial lattice spacing and 



SL = r_ 2 S >lX ' (43) 



P2 1 n - p£ 

2m 



< fl. 



In computing the integral over the spatial momentum on the lattice one usually assumes 
that the cutoff is far from the internal scale of the model, takes the limit a s — > and 
expands the lattice propagators. We regain the continuum propagator in this manner and 
the higher order terms in this expansion, the irrelevant, higher order derivative operators are 
formally suppressed by a positive power of the lattice spacing. One is tempted to ignore them 
altogether. But this is not always allowed because the ultraviolet divergences of the loop 
integration might generate so strong singularities that quantities which seemed to disappear 
turn out to be finite or even divergent as a s — > 0. 

One can see that no problem arises if one considers observables in a renormalizable 
model without anomaly. We can always carry out the substractions in the loop integrands of 
a renormalizable theory in such a manner that the resulting integrals are finite. An anomaly 
appears when "by accident" a graph with non-negative primitive degree of divergence given 
by the power counting happens to be finite and receives no subtraction during the renor- 
malization process. Thus in a renormalizable, non-anomalous model all loop integrals of 
the renormalized perturbation expansion are finite and have negative primitive degree of 
divergence. These integrals converge uniformly as the cutoff is removed and the order of 
the integration and the limit a s — > can be interchanged. By setting a s = in the inte- 
grands one eliminates all lattice artifacts and the loop integrals which follow reproduce the 
continuum perturbation expansion. 

Returning to our non-relativistic model, the lesson of the argument about the suppression 
of the lattice artifact is that the continuum propagators can safely be used for quantities 
which stay finite as a s — > 0. This is enough to simplify the computation of several important 
quantities such as the imaginary part of the self energies. For other observables, the cutoff 
a s must be kept finite and the ultraviolet details of the model remain important. 



VI. ELECTRON PROPAGATION 

In this section we calculate the first two order contributions to the electron self energy 
in two dimensions. The goal of this calculation is to show that the diagrams with at most 
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two impurity lines atached to each vertex agree in the present approach and those of the 
Edwards model can be brought into equivalence after the rescaling of the impurity density. 

We consider the case A = for simplicity when the Feynman rules can be read off from 
the generating functional 



xg i J dtdx[J* (f>+(j>* J+jti/>a+iptja.] e -iji J dtdxip' a ijj a cj)* <j> 

, -Hf Mir J 5 S 6 

det[(iG a ) _1 ] det[(iGo) -1 ]e R J «t (M) 

xe « J dtdxdt 1 'dx'[J*(x,t)G a (x,t,x' \t')J(x' ' ,t')+'p (x,t)G (x,t,x' \t')j(x' ,t')] 



The electron self energy in the leading order is computed in Appendix [E 1| and Eq. ( |E5| ) 
yields 

= pg. (45) 
The two 0{g 2 ) contributions, shown in Figs. [I], § and computed in Appendix E2| give 



and 



ImZ?\uj,p) = --g 2 p(l - p)no(jJL), (46) 



ImZf\uj,p) = -g 2 p(l -p)j {29{u) - 9(co + p)) , (47) 



respectively where e = + , and 



no(p) = ^- (48) 



There are peculiarities about the graph in Fig. 0. First, it involves the feedback of the 
electron dynamics on the impurities and it is thereby suppressed in the replica method. What 
is the role of this contribution in our model? The other peculiarity is that the meaningless 
1/e in Eq. (fIBD indicates an infrared divergence, expected to appear in zero dimensional 
systems. We show that the solution of the latter problem is the answer for the question 
raised before, namely it removes this particular graph altogether from our model. The point 
is that the divergence can be removed by the usual ring diagram, RPA resummation, by 
taking into account the "screening" of the electron density in the impurity propagator so 
that 



H a (t) = a' J dx(f) j (x,t)(j)(x,t) -> H a+2 Tgn (ti)(t) = (a' + gn (p)) j dx<j) j (x,t)(j)(x,t). (49) 
This procedure yields 

iG a+2Tgno{ll) {x^x\t') = 8{x - x i) e -^T 9 n ( P )Kt-t') (50) 

X \6(t - t' - T]) (l - n a+ 2T fl no(At)) ~ ^ ~ t + r l) n a+2Tgn (n)} ■ 

The correction to the electron propagator is therefore 
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1 Z" 71 " Cj f 

— J dae ian detKiGa^TgnoM)' 1 }!-^ J dxiGo{x 1 ,x)iG (x,x 2 )iG a+ 2 Tgno ( fJ ,) 

r -, TV— 1 

g P dae ian ( 1 + e -*O+2Tgn (/i))) e -i(a+2Tgn (^)) 



ha? s dae iaN [1 + e -^ a+2T 9 n oW) 

dxG (xi, x)G (x, x 2 ) 



2V-1 

e -i(a+-ZTgn (n)) f 

dxG (xi,x)Go(x, x 2 ) 



N 



g \n-l 



ha? s fN 
V n 

■ 99 



dxG (xi,x)G (x,x 2 ). (51) 



The result of such a partial resummation is finite and agrees with (|4~5|) . This is because G a 
is a periodic function of a and the integration, carried out over the length of a period, is 



invariant under the shift a — > a + c. Thus the only effect of the renormalization (49) is the 
omission of the contribution of the graph of Fig. |T] to the RPA ring diagrams altogether. 
In other words, the discretness of the electron number removes the feedback of the electron 
dynamics on the impurities. 

The complete 0(g 2 ) self energy after resummation is 

Im^ 2 \u) = -g 2 p(l - (29(co) - 9(cu + //)) , (52) 

which is independent of the momentum p because the impurity dynamics is ultralocal, the 
pointlike impurities do not propagate. The imaginary part of the self energy is related to 
the mean free path, in particular for the electrons above of the Fermi level, w>0we have 

77? Tl 

Im£H = -g 2 p(l - p)- + 0(g 3 ) « --, (53) 

2 2r 

where r is the relaxation time. Recall that the dimensionless density p G [0, 1] and 
JmS = O(p), the electron relaxation time diverges at low impurity densities. The anal- 



ogous expression of the Edwards model can be obtained from (]52|) by the replacement 



p(l -p)^p. (54) 

Our expression reflects the fact that the electrons do not scatter on the completely filled up 
impurity system (Pauli blocking), a phenomenon which is neglected when the effects of the 
impurities is represented by a static potential only. 

Let us now consider the propagation of a particle- hole pair with energy- momenta u + ,p + , 
respectively, and write the corresponding amplitude as 

G(uj + ,p + ,uj_,p_,uj' + , p' + , wl, pL) = iG (w + , p+)iG {u)' + , p' + )iG Q (w_ , p- )iG (w'_ , pL) 

xS(p + — p\ — p- + p'_)E. (55) 

In the leading order, before the a integration we find 

E a = (ig) 2 X! iG a {<jj)iG a {ijj')5{uj + — oj' + + oj — u)')5(u' — uj + oj'_oj_) 

= -(ig) 2 2n(l-p a )p a . (56) 
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The a integration yields 



(^) 2 2tt(1 - p)p, 



(57) 



a result which agrees with the Edwards model after the replacement (|54|). The insertion of 
i impurity bubbles at different sites and the a integration yield the factor 



Since the particle-hole lines are the usual ones in the rest of the diagram, this leads to the 
equivalence of the resummation over the ladder or the maximally crossed graphs, as well. 
In the light of the results of ref. || this indicates that the localisation-delocalisation phase 
structure is similar in these descriptions. 

Not all diagrams of the local, Gaussian potential model (HD can be brought into equiv- 
alence with our model. There one finds contributions with more than two impurity lines 
at a coordinate space location. This can not happen in our model where the fermionic 
statistic of the impurity cancels these diagrams. But it is easy to see that the remaining 
diagrams, where only two impurity lines are attached at each spatial lattice site agree after 
the change (|54l) . Furthermore, the disconnected contributions simplify for these graphs in 
the thermodynamic limit. The point is that the impurity propagators factorise and the 
disconnected parts drop out as in the usual, annealed averages. The resulting impurity 
propagators at different lattice sites give the same contributions when the thermodynamic 
limit is performed. This result is actually expected if the canonical and the grand canonical 
ensemble are equivalent. 

There are furthermore diagrams in our model which have no counterpart in (|lj). These 
graph contain internal fermion loops and some of them cancel after resummation. But there 
are remaining non-vanishing contributions, e.g. the correlation of the potential is ultralocal 
and unchanged by the electrons for ([!]) but the electron dynamics generates a non-trivial 
correlation function for the impurity density (0*0) in our model, the remnant of the annealed 
integration over the field 0. Such an electron dynamics generated correlation is absent in 
the Green functions for the impurities, the correlation functions of the field variable 0* or 
remain local due to the gauge symmetry ([|). It remains to be seen if the feedback of the 
electron dynamics generated by the resummation of the irreducible vertices obtained in the 
higher loop aproximation cancels or not. 



The quenched averaging over static impurities is achieved within the second quantized 
formalism by means of an annealed averaging procedure where the impurity motion is slowed 
down by setting the impurity mass equal to infinity. This description corresponds to a system 
of electrons and static impurities in equilibrium. Despite the feedback of the electrons on the 
impurity dynamics the localisation-delocalisation phase structure in the present approach is 
expected to be similar to the phase structure of the Edwards model. 



W - p)\ 



(58) 




VII. CONCLUSIONS 
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The problem of the usual methods for disordered systems is their inability to deal with 
strong interactions. However our model is well suited to non-perturbative methods. The nu- 
merical simulation in lattice regularization is feasible since one can easily construct stochas- 
tic sampling algorithms in our canonical ensemble at finite temperature. Another natural 
non-perturbative method available is the functional renormalization group approach in the 
internal space [III] which can treat the finite density fermionic systems in a simple manner. 
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APPENDIX A: LATTICE REGULARIZATION 

The regulator applied in this work is a lattice in space-time with spatial and temporal 
lattice spacing a s , and a t , respectively. The action is written as 

S = lJ GMrL^a;,*),^^ (Al) 

with S = Si + S e + S', 

h J dtdx(f)\x,t){idt — a)(j){x,t) — > Si = 0t f(idt — 

£,t 

J dtdx^(x,t)(id t + + n)il){x,t) ^S e = Y,yl,0t + + /*)4,f, 

x,i 

-gh f dtdx^{x,t)^{x,t)<p\x,t)<P{x,t) ^S' = -gY,^{AA,Ah ( A2 ) 



a)0 £ t, 



x,t 



where x- 7 = n^a s , j = 1, ■ • • , d, t = £a t and the following lattice quantities were introduced 
in 2D (fizf = a s (j)(x,t), (the field variables with space-time coordinates in the subscript are 
dimensionless) a = a t a, rh = a s m, h = a t h/a s , ip^f = a s ip(x,i), g = g/a s , and ft = a t fi. The 
finite difference operators are defined as d T f n/ = f n/ - / n ^_i and Af n>e = Yh=i fn+e z ,e + 

fn—ei,l 2f n £. 
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APPENDIX B: OPERATOR FORMALISM 

The dynamics generated by the free impurity hamiltonian 

X 

where the static impurity density is denoted by 

is simplest to discuss in the operator formalism. The impurity propagator is written 
iG a (x,t,x',t') = - 2 C n .,.. l , l ,r = iE(»l r [ e "* J ^* ff "' ( * ) te*)# t (^0]|n>, 

a s * {„} 

where 

in>=n (&)*»■ 

X 

and M is a normalization constant. The Hamiltonian ( pT|) gives 

w 

M 

Since the lattice sites are decoupled we find immediately 

E { „}(n|e- iQ '( T -')4e- ia '(*-*')4e--'(''+ T )|n)= e~ ia ' ^5^,(1 + e~ ia ) N ~\ 
E { n}(n\e- ia ^ T - t '^ie- ia '^'- t ^ £ e- M ^\n)= e'™'^ 8^,(1 + e- ia ) N ~ l e~ iaT , 

where N denotes the number of sites. The normalization constant is 

M = Aet[iG a }- 1 = (l + e~ ia ) N , 

yielding 

'('„.,,,:r = e- iQ '('-''^(l + e^Y 1 (6(t - t') - 6(t' - t)e~* a ) 
on the lattice or 

iG a {x, t, x', t') = 5(x - x')e- ia ' {t - tr) (6(t - 0(1 - n a ) - 6{t' - t)e~ ia n a ) 
in the continuum where 

1 

n a = - — . 

1 + e ia 
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APPENDIX C: FUNCTIONAL APPROACH 



The interactions are easier to take into account in the path integral representation. To 
make sure that all singular features of the limit of diverging impurity mass are properly 
kept we first quickly reproduce the free impurity propagator in this representation. The 
amplitude 



Tre J - T 



i C T dtH a ,(t) 



(CI) 



can be written as a functional integral over Grassmannian coherent state configurations 
which are antiperiodic in time by inserting the resolution of the identity 



J] d(f)*(x, t)d(j){x, t)e~ I «***•(*.*)*(».*) t) (<j)(x, t) | , 



(C2) 



into ( |Cip at different times t m , 



Z 0a = Jim 



N t 

n n< 

m=— N t +l x 



[X 



N, 



x II (0(x,t m )|e- ia '^'«|0(x,t m )) 

m=-Nt+l 

lim / TT n^^^^^e"^ 1 ^'^^ 

Nt-ooJ m= _ Nt+1 x 

•SmL-ATt+a/ dx <t>*(. x t™)(l-iatot')cl>(x,t m -i) - J dx<f>* (x,t- Nt+1 )(l-ia t a')<p(x,t Nt ) 
N t 



xe 



IT [ lim / IT d<f>*(x,t m )d<f>(x,t m ) 



m=-N t +l 



-,N t 



g Z— /m,n— — Nf 



+1 <f>*(x,t m )S mn <f>(x,t n )-\ 



(C3) 



where the matrix S is 



/I o 

-A 1 ■•• 

-A 1 



S 



A \ 
... o 
... o 



(C4) 



with A = 1 — iatot 



\ -A 1 / 

1 — ia/2N t . The Grassmann integration gives 



lim det S 

N t ->oo 



x 

The propagator is given for t m > t m i by 



TT lim 

J - J - iV t ->oo 

X 



1 + A 



2N t 



1 + e 



N 



(C5) 
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iG a % x'(t m , t m i) — S x x i lim (S ) 

N— >oo \ / mm' 



S T r' lim 



TV^oo 

j^m—m' 



N ^oo 1 _|_ A 2Nt 

S^e-^-^fl-^—), (C6) 



and for i m < £ m / by 

zGq,* x'(t m , t m i) = 8x,x' lim [S ) 

/V+— >oo \ / mm' 



(5-?. lim 



_^2iV t +TO-m' 



A^oo 1 + A 2iV ' 

-^e- fa * (C7) 

1 + e za 



in agreement with (|E8|). 

The Fourier transform of the propagator in time is 

G a (x, x', to) = f d(t - t'y^-^Gaix, t, x', t') 



8(x - x') lim ( 1 Ha + - ) (C8) 



in the continuum and 



/I — 77 71 \ 

GoA* W = lim e"-" — ry^- + _ (C9) 

m^o+ y ^ - a' + ie ™ - a' - ie j 

on the lattice where the infinitesimal quantity r\ was introduced to take care about the point 
splitting. 

We shall need later 

00 rl(,i 



— iG a (x, x; uj)e 

-00 2,71 



r)^0+ J — 00 2tt 

i lim r°° ( 1 ~ n a ^ 



a? s v^o + J -00 2ir \cu — a' + ie uo — a' — ie 



i . f°° du n n e iuJV 
— hm / — 



a? s v^o + J -00 2n lu — a' — ie 
1 

a 



2 (l + e- ia )- l e~ ia (CIO) 



which, after integration over a leads to 



1 r dae ina (1 + e - "*^ 1 e 
0,2,0+)- J - 1 + ' 



aj J^dae ina (l + e- ia y" 
n 

p. (Gil) 



18 



APPENDIX D: EFFECTIVE THEORY 



Consider now the transition amplitudes of the impurity-electron coupled system in the 
path integral representation 

(*/,e,*/ I i|e- i ^|*i > e,*i,i> = II / #^,«#*,f,a#iA,f 

x,t 

of the model (H) on a discretised space-time lattice, where 



ft 



2m 
A 



x+a s i,t-a t , 



a ^Px—a a ii—a t ,a ^^Px,i—a t ,a) 



(D2) 



~ 4>x,t-a t ) + 9^ii,a^x,i-a t ,a0l,i4>x,i- 



The wave functional of the initial and the final states are ^i >e [ip*, ip], ^$[0*, fa, ^% e [^*, 'fa, 
and 

Notice that the impurity field is coupled in the time direction only and the functional 
integral over fa and <fi decouples into a single site, quantum mechanical problem. This 
suggests the introduction of an effective theory obtained after integration over the impurity 
field, 



n 



x,t 



This is the path integral of a two level system with hamiltonian H t = J t , where 

J t = 9^ti,a^x,t-a t ,a- 

The matrix elements 

(0|e-^|0) = 1, 
(l\e' l ^ H \l) = T 



(D3) 



(D4) 



lie 

L t 



-iJt 



agree with the path integral 

(n|e-^|n) = (0^rn 



(D5) 



(D6) 



where fa and 0/ stand for the first and the last integration variable in time, the arguments 
of the initial and final single site wave function 
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(D7) 



In fact, (P6|) can be checked by either direct computation of the fermion determinant or by 
expanding the integrand in the quasi-local term (f)*(j) t _i. One finds the effective action 



for the states 



2m ^ Y *>w Yl 



x+a s i,t—at 



,a + ^x 



x—a s i,t—at 



n ^^x.i—at.a) 



x,t—at,a 



Hi 



n 



'x i 



(D8) 



(D9) 



which belong to the occupation number configuration {n x }, = 0,1. Notice that the 
impurity induced interaction is as long range in time as possible, i.e. time independent, and 
ultralocal in space due to the conservation laws originating from the symmetry (H). 

The time ordered expectation value of an observable O of the electrons can easily be 
obtained as a summation over the impurity occupational number configurations n$, J2x n x = 



where 



«*/,e|T[0]|^, e » = 7 EII d^J^a 

P -K} £,t 

** e [^*,V']* i)e [Vi*^] e -i^/^^C», 



x,t 



(D10) 



(Dll) 



APPENDIX E: ELECTRON SELF ENERGY 



1. O(g) 



The a-dependent electron propagator can be written in terms of the free electron and 
impurity propagators Go and G a , respectively as 

(^(x 1 )^(x 2 )) a = det[{iG a )- 1 } det^Go)- 1 ] 

x< iGo(xx, x 2 ) I — ijr J dxiG a (x,x)iGo(x,x) 



■ 9 



n 



dx / dy(iGo(x,x)iGQ(y,y)iG a (x,x)iG a (y,y) 
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-iG a (x, x)iG a (y, y)iG (x, y)iG (y, x) 
-iG a (x, y)iG a (y, x)iG (x, x)iG (y, y)) 

J dxiGo(xi, x)iGo(x, x 2 )iG a (x, x) 



x 



/ dyiG (y,y)iG a (y,y) 

■ 9\' 



+ [ i t) J dX J d y{ iG o( x ^ x ) iG o(x,x 2 )iG a (x,y)iG a (y,x)iG (y,y) 
-iG (x 1 ,x)iG (x, y)iG (y, x 2 )iG a {x, y)iG a {y, x)) | + 0(g 3 ), (El) 

where the compact notation (x, t) — > x was introduced. For the normalization we need 

Z a = J D[^}D[4>]D[(f)*}D[(f)}ei J ^xdt[^{i%d t -na')4>+^{ind t +^+m 

x ^1 - i| J dtdxip\x, t)if>(x, t)(f)*(x, t)(f)(x, t) + 0{g 2 )j 

= det[iG o r 1 (l + e- M ) 

x|l-i|y dxdtiG (xt;xt)iG a (xt;xt) + (% 2 )j. (E2) 

The complete propagator 



starts as iGo(xi, x 2 ) in The non-trivial O(g) piece in the numerator is 

dae tan det[(iCr a ) -1 ]i^ / dxiGo(xi, x)iGo(x, x 2 )iG a (x, x) 



x 



X ~ i \\ d y iG o(y,y)iG a (y,y) 



= % 



■ 9 
Tial 



• 9 r vn i /iV-2' 



Taking into account the normalization we find 

(ip{xi)i) ] {x 2 )) =iG (x 1 ;x 2 ) 



+i ~ha? ' dxG (x 1 ,x)G (x,x 2 ) 



N-l\ . a „ r .„ , (N -2 



n- 1 



= i(?o(^i; ^2) + y dxG (xi; x)G (a;; x 2 ) + 0(g 2 



(E3) 



(^-2) J dxG o{xi,x)G Q {x,x 2 ). (E4) 



+ 0(g 2 ) 



(E5) 



in the thermodynamic limit, (N,n — > 00, ra/iV — > p). The comparison with the Schwinger- 
Dyson equation equation gives the self energy 
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£« = pg, (E6) 

a simple shift in the energy. In general, the diagram with I impurity bubbles produces the 
factor 

(N-t\ 

\n-£j N\(n-£)\ £ 
(N\ n\(N-£)\ ' ' 
V n 



2. 0(g 2 ) 

The graph shown in Fig.|I] is 

I a i = J dxdtdydt'iGo(xi, t\\ x, t)iGo(x, t; X2, t 2 )iGa(y : t'; y, t')iG a (x, t; y, t')iG a (y, t'; x, t) 

= — J dx J dt J dt'iG (xi,ti;x,t)iG (x,t;x2,t2)iG (y,t';y,t') 

x [0(t - f) (1 - n a ) - 6(t' - t)n a ] [9(1/ - t) (1 - n«) - 0(t - t')n a ] 
i 5 r diut f duj 2 f duj 3 f d 2 p! f d 2 p 2 e^i-^-i^h-t^ 



a 2 J 2tt J 2vr J 2vr J (2n) 2 J (2vr) 2 fa _ |i + ^ + 



2m 
2 



x ^ e ***u(2_J^ + _^) . (E8) 

The integration over CJ3 is divergent 



27T \Cg>3 + Z6 Co> 3 — 26, 

where e — > + . The projection onto the subspace with n impurities is obtained by the 
integration 

Z n e ian det[{iG a y l ]{l - n a )n a _ e^-^l + e ~ ia ) N - 2 
detf^)- 1 ] ~ /^^^(l + e-^ 

'iV-2' 



n 

p(l-p). (E10) 



The integration over u 2 and j?2 gives 



d 2 p 2 f duj 2 1 /" d 2 p 2 .. f duj 2 e" 



lim 



iu>2t 



(2n) 2 J 2tt ^ _ + ^ + i<% 7 (2tt)W7 2tt ^ _ gL + ^ + i% 

= J (2^^-2^ 

= zn (/i). (Ell) 
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The corresponding contribution to the self-energy, obtained after removing the external legs 
is 

4 2 V,p) = - Vp(l - p)no(ji). (E12) 
The graph of Fig. |2| gives 
I a 2 = J dxdtdydt'iGo(xi,ti; x, 1)iGo(x, 1; y, t')iGo(y, t'; X2, t2)iG a (x, t; y, t')iG a (y, t'; x, t) 

= — J dx J dt J dt'iGo(xi,ti; x,t)iGo(x,t; x,t')iGo(x,t'; X2,t2) 

x [9(1 - t') (1 - n a ) - 9(t' - t)n a ] [9(1? -t)(l- n a ) - 9(1 - l')n a ) 
i r duox r duo 2 fdu 3 f d 2 p x r d 2 p 2 e *i{*i-**)-*>*fr-**) 



a 2 J 2tt J 2vr J 2vr J (2n) 2 J (27r) 2 ( ^_|l + yU + ^ i)2 
1 / 1 — n n n n 



x- 



( 1 — n„ n n \ . 

x - -+ ■ E13 

\ou 3 + u) 2 - UJ\ + ie uj 3 + u 2 - 0J1 - xej 

The integration over uo 3 of the terms between the brackets yields 

z(l - n a )n a ( — — -) = 2tt(1 - n^nj^x - u 2 ), (E14) 

V uj\ — LU2 + te io\ — LU2 — lej 



thus 



i fdut r d 2 p 1 r d 2 p 2 e ipi{xi - X2) - iuJl{tl - t2) 1 

Ia2 = / TT" / tttt? / 77TT? 7a 12 1 1 - ™a)™o 



2W (2^ (2tt)* (a7i _ £ + ^ + i(%)2 ^ _ Sl + ^ + i<% 

(E15) 

The integration over p 2 gives 

/ ggg 1 r ggg (pp 1 f s 1 fE16) 



where P.P. stands for the principal part and produces only a shift in the energy. This shift, 
being divergent as a s — > is not computable without specifying the details of the ultraviolet 
cutoff. But the more important imaginary part is finite and can be decomposed in the 
following way 



— 1 



f d 2 p 2 = i— ( F - deit5(u + -) 

J (2ir) 2 ( 0Jl _£ + fi) 2 + 5 2 *2vrU J„ J " } 

771 

= i- {[9(00, + //) - 9(u x )\ - . (E17) 

The term containing 9(ui + /x) — 9(uo\) represents the contribution of the Fermi sphere. The 
last term stands for the contributions of the excitations above the Fermi surface. Their sum 
is proportional to the free electron density of states, 
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N (E) = —2 (20(25 -p)- 9(E)) (E18) 

Z7TO 

which is positive when E is above the Fermi level and negative for < E < p. The final 
result 

t ' Tfl 

imi 2 (u, P ) = - 2 p(i - p)- — — rW"9 W w + ^ - 2 ^ • ( E19 ) 

gives the imaginary part of the self-energy 

ImY$\u,p) = -ig 2 j P (l - p) {26{u) - 6{uo + p)) . (E20) 



24 



FIGURES 




i t 

\ / 

> ^ > 

FIG. 1. Energy independent 0(g 2 ) contribution to the electron self energy. 




FIG. 2. Energy dependent 0(g 2 ) contribution to the electron self energy. 
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